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Abstract 

We prove optimality conditions for different variational functionals 
containing left and right Caputo fractional derivatives. A sufficient condi- 
tion of minimization under an appropriate convexity assumption is given. 
£> An Euler-Lagrange equation for functionals where the lower and upper 

£ | bounds of the integral are distinct of the bounds of the Caputo derivative 

CO ■ is also proved. Then, the fractional isoperimetric problem is formulated 

CTN ' with an integral constraint also containing Caputo derivatives. Normal 

and abnormal extremals are considered. 
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1 Introduction 

5-H 



c3 



Fractional Calculus is the branch of Mathematics that generalizes the deriva- 
tive and the integral of a function to a non-integer order. The subject is not 
recent and it is as old as the calculus itself. In a letter dated 30th September 
1695, Leibniz proposed the following problem to L'Hopital: "Can the meaning 
of derivatives with integer order be generalized to derivatives with non-integer 
orders?" Since then, several mathematicians studied this question, among them 
Liouville, Riemann, Weyl and Letnikov. There are many fields of applications 
where we can use the fractional calculus, like viscoelasticity, electrochemistry, 
diffusion processes, control theory, heat conduction, electricity, mechanics, chaos 
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and fractals (see some references at the end, e.g., [rUlTT71l?Tll331IMll371l5Ull^ ) . To 
solve fractional differential equations, there exist several methods: Laplace and 
Fourier transforms, truncated Taylor series, numerical methods, etc. (see [5] and 
references therein) . Recently, a lot of attention has been put on the fractional 

calculus of variations (see, e.g., mi3HigiS[Iim[H|5Dl|2Sl|551|251|25]. W& ah ° 
mention [7 j , were necessary and sufficient conditions of optimality for functionals 
containing fractional integrals and fractional derivatives are presented. For re- 
sults on fractional optimal control see [2"IIT5] . In the present paper we work with 
the Caputo fractional derivative. For problems of calculus of variations with 
boundary conditions, Caputo's derivative seems to be more natural, since for a 
given function y to have continuous Riemann-Liouville fractional derivative on a 
closed interval [a, b], the function must satisfy the conditions y(a) = y(b) = [8]. 
We also mention that, if y(a) — 0, then the left Riemann-Liouville derivative is 
equal to the left Caputo derivative. 

The paper is organized as follows. In Section [5] we present the necessary 
definitions and some necessary facts about fractional calculus. Section [3] is 
dedicated to our original results. We study fractional Euler-Lagrange equations 
and the fractional isoperimetric problem within Caputo's fractional derivative 
context for different kinds of functionals. We also give sufficient conditions of 
optimality for fractional variational problems. 

2 Preliminaries 

2.1 Review on fractional calculus 

There are several definitions of fractional derivatives and fractional integrals, like 
Riemann-Liouville, Caputo, Riesz, RieszCaputo, Weyl, Grunwald-Letnikov, 
Hadamard, Chen, etc. We will present the definitions of the first two of them. 
Except otherwise stated, proofs of results may be found in |18| . 

Let / : [a, b] — >• R be a function, a a positive real number, n the integer 
satisfying n — 1 < a < n, and T the Euler gamma function. Then, 

1. the left and right Riemann-Liouville fractional integrals of order a are 
defined by 

alSttx) = -^ [ X (X - t) a - l f(t)dt, 



and 



respectively; 



T(a) 
1 



xl?f(x) = t^T / (* - X) a - l f{t)dt, 



2. the left and right Riemann-Liouville fractional derivatives of order a are 
defined by 



and 



dn rn- at ,„, (-1)" d n 



x D?f(x) = (-1)™^— x I"- a /(a;) = =i — '—--r- / (t - xY 1 - - 1 f(t)dt, 

respectively; 
3. the left and right Caputo fractional derivatives of order a arc defined by 

An 1 j-x 

c D a x f(x) = a i:- a —f{x) = T{n _ a) J a (x-tr-^f^(t)dt, 

and 

An 1 fb 



respectively 

There exists a relation between the Ricmann-Liouville and the Caputo frac- 
tional derivatives: 

c a D«f(x) = a D»f(x) J2 r / fc)(a L (* - «) fe_a 
and 

to r ( fc - « + 1) 
Therefore, 

if f(a) = f'(a) = ... = /(""^(a) = 0, then c a D a J{x) = a D a J{x) 
and 

if f(b) = f'(b) = ... = .^-^(b) = 0, then °D?f(x) = x D%f(x). 
These fractional operators are linear, i.e., 

V(ii.f(x) + ug{x)) = nVf(x) + vVg(x) 1 

where V is a D" , X D^,^D X , x D^, a I x or X I", and /u and v are real numbers. 

If / G C n [a, 6] , then the left and right Caputo derivatives are continuous on 
[a, b]. The main advantage of Caputo's approach is that the initial conditions 
for fractional differential equations with Caputo derivatives take on the same 
form as for integer-order differential equations. 

Some properties valid for integer differentiation and integer integration re- 
main valid for fractional differentiation and fractional integration; namely the 
Caputo fractional derivatives and the Riemann-Liouville fractional integrals are 
inverse operations: 



1. If / € L OQ (a, b) or / € C[a, 6], and if a > 0, then 

c a D a xa I^f{x) = f{x) and °D? w I?f(x) = f(x). 



2. If / e C">, 6] and if a > 0, then 



/ (fc) («) 



ai: c a D a j{ X ) = f( X ) J2 hr 1 ^ fl ) 



fe=0 



fc! 



and 



fc=0 



fc! 



We also need for our purposes integration by parts formulas. For a > 0, we 
have (cf. @]) 



b i-b 



g(x) ■ ^D«f(x)dx = / f(x) ■ x D%g{x)dx 



J2[*DZ +j - n g(x)- x DZ- l -if(x) 



3=0 



and 



b r b 



j(x) ■ ^DZf(x)dx = / f{x) ■ a D a x g{x)dx 

J a 
n-1 

+ E [(-l)" +3 «^ +r "s(^) • a D:~ l - j f(x)] b a , 

where a D%g{x) = J x k gix) and x D^g(x) = x I^ k g(x) if fc < 0. 
Therefore, if < a < 1, we obtain 

/ 0(aO • C a D a J(x)dx = f f{x) ■ x D«gix)dx + [J^gix) ■ fix)} \ (1) 

J a J a 

and 

f g(x) ■ °;DUix)dx = I fix) ■ a D«g(x)dx - [j^g (x) ■ fix)} \ . (2) 

Moreover, if / is a function such that f(a) = f(b) = 0, we have simpler formulas: 

/ g(x) ■ C a D«f{x)dx = f fix) ■ x D%gix)dx (3) 

J a J a 



and 



b pb 



j(x) ■ ^D^fix)dx = / f(x) ■ a D a x gix)dx 



(1) 



Remark 1. As a goes to 1, expressions ([IJ and ^ reduce to the classical 
integration by parts formulas: 

g(x) ■ f'(x)dx = - I f(x) ■ g'(x)dx + [g(x) ■ f(x)] b a 
and 

r b rb 

9{z) ■ {-f'{x))dx = / f{x) ■ g'(x)dx - [g(x) ■ f{x)] b a , 

J a 

respectively. 

Remark 2. Observe that the left member of equations ([3]) and ((H) contains a 
Caputo fractional derivative, while the other one contains a Riemann-Liouville 
fractional derivative. However, since there exists a relation between the two 
derivatives, we could present this formula with only one fractional derivative, 
although in this case the resulting equation will contain some extra terms. 

2.2 Fractional Euler— Lagrange equations 

From now on we fix < a, ft < 1. Also, to simplify, we denote 

[y]{x)~{x,y{x), c a D^y{x)^Dly{x)). 
Let 
"E^ — < y : [a, b] — > R | ^^xV an d x^bV exis t and are continuous on [a, b] > . 

Definition 2.1. The space of variations Var(a, b) for the Caputo derivatives 
is defined by 

c Var(a, 6) = {/i£ %E% \ h{a) = h(b) = o} . 

We now present first order necessary conditions of optimality for functionals, 
defined on "E^ , of the type 

J(y) = f L[y](x)dx. (5) 

J a 

We assume that the map (x,y,u,v) — > L(x,y 7 u,v) is a function of class C 1 . 
Denoting by c^L the partial derivative of L with respect to the ith variable, 
i = 1,2,3,4, we also assume that d^L has continuous right Riemann-Liouville 
fractional derivative of order a and d^L has continuous left Riemann-Liouville 
fractional derivative of order ft. 

Definition 2.2. We say that y is a local minimizer (respectively local max- 
imizer) of J if there exists a 5 > such that J(y) < J(yi) (respectively 
J(y) > J{yi)) f° r au J/i such that \\y — yi\\ < 5. 



In [T] Agrawal considers the problem of finding extremals for functionals 
containing left and right Riemann-Liouville fractional derivatives of the form 

,& 
J(y) = / L(x,y{x), a D%y(x), x D p b y{x))dx 

•J a 

and he derived an Euler-Lagrange equation for an extremum y of J, subject to 
the boundary conditions y(a) — y a , y(b) = y\,\ 

d 2 L + x D£d 3 L + a D^diL = 0, for all x e [a, b]. 

In [9] a new type of functional is studied, in case where the lower bound of 
the integral do not coincide with the lower bound of the fractional derivative: 

J*(y)= L(x,y(x), a D%y(x))dx, 



where [A, B] C [a,b]. We also mention |3], where an Euler-Lagrange equa- 
tion and a transversality condition are given, for functionals with left Caputo 
derivatives and a boundary condition on the initial point x = a. 

Theorem 2.3 ( [3]). Let J be the functional defined by 

rb 



J(y)= / L(x,y(x)^D^y(x))dx. 

J a 

Let y be a local minimizer of J satisfying the boundary condition y(a) = y^. 
Then, y satisfies the following conditions: 

2 L + X D£3 3 L = (6) 

and 



xJ-b 



"9 S L\ b = 0. (7) 



Note that the Euler-Lagrange equation ([5]) contains a Riemann-Liouville 
fractional derivative, although J has only Caputo's derivative. Also, the transver- 
sality condition ([?])■ in general, contains fractional derivative terms. Thus, in 
order to solve a fractional variational problem, it may be required fractional 
boundary conditions. This result is then proven for functionals with higher 
order fractional derivatives. In [5] functionals with the left and right Caputo 
fractional derivatives are considered. We include here a short proof. For more 
on the subject we refer the reader to [TTinBll^B] . 

Theorem 2.4 ( [5]). Let J be the functional as in ^ and y a local minimizer of 
J satisfying the boundary conditions y{a) — y a and y(b) = yb- Then, y satisfies 
the Euler-Lagrange equation 

d 2 L + x D«d 3 L + a D%L = 0. (8) 



Proof. Let e be a small real parameter and 77 6 c Var(a, b). Consider a variation 
of y; say y + erj. Since the Caputo derivative operators are linear, we have 

,b 
J(V + «/) = / L(x, y + e V , c a D a x y + e^D^t], c x D>ty + e^D§r,)dx. 

We can regard J as a function of one variable, J(e) = J(y + eij). Since y 
is the local minimizer, J attains an extremum at e = 0. Differentiating J(e) at 
zero, it follows that 

f b r i 

/ d 2 L V + d 3 L^D^r] + d 4 L^D^ri dx = 0. 

J a 

Integrating by parts, and since 77(a) = 77(6) = 0, one finds that 

,6 

/ [d 2 L + x D£d 3 L + a D%L] r,dx = 

J a 

for all 77 G c Var(a, b). By the arbitrariness of 77 and by the fundamental lemma 
of the calculus of variations (see, e.g., [32j p. 32]), it follows that 

8 2 L + x D£d 3 L + aDldiL = 0. 

□ 

When the term x D^y is not present in the function L, then equation ([8]) 
reduces to a simpler one: 

d 2 L + X D«8 3 L = 0. 

Moreover, if we allow a = 1, and since in that case the right Riemann-Liouville 
fractional derivative is equal to —d/dx, we obtain the classical Euler-Lagrangc 
equation: 

d 2 L - ^-d 3 L = 0. 
dx 

3 Main results 

Our first contribution is to generalize Theorem 12.41 (cf. Theorem 13.11 below). 
We consider a new type of functional, where the lower bound of the integral is 
greater than the lower bound of the Caputo's derivative, and the upper bound 
of the integral is less than the upper bound of the Caputo's derivative. Because 
of this, we can not apply directly the integration by parts formula and some 
technical auxiliary procedures are required. A similar problem is addressed for 
the Riemann-Liouville fractional derivatives in [9]. 



3.1 The Euler-Lagrange equation 

We consider the functional 



J*(y) 



CnP„ 



L(x,y(x)^D«y(x)^D£y(x))dx 



(9) 



where (x, y, u, v) — > L(x, y, u, v) e C 1 and [A, B] C [a, b]. 

Theorem 3.1. If y is a local minimizer of J* given by (0J), satisfying the 
boundary conditions y{a) = y a and y(b) = yb, then y satisfies the system 

d 2 L + x D a B d 3 L + A D^d A L = for all x € [A, B] 
x D B d 3 L - x D a A d 3 L = for all x e [a, A] 

A D^.d 4 L - B D^d 4 L = for all x e [B, b] 

Proof. Let y be a minimizer and let y = y+erj be a variation of y, r\ £ Var(a, b), 
such that n{A) = n(B) = 0. Define the new function J*(e) = J*(y + en). By 
hypothesis, y is a local extremum of J* and so J* has a local extremum at e = 0. 
Therefore, the following holds: 



= 



d 2 Ln + d 3 L^D^ + d 4 L^D p b r, 

B 



d 2 L rjdx 



dx 

A 



f d 3 L^D^r,dx- f d 3 L c a D a x nd 

J a J a 



+ 



fb pb 

I d 4 L^D^ndx- / diL^D^ndx 

J A ' JB 



I A JB 

Integrating by parts the four last terms gives: 







B 



d 2 L ndx 



B 



n x D B d 3 Ldx 



n x D A d 3 Ldx 



pb pb 

/ n AD^diLdx — I n sD^diLdx 

J A JB 



d 2 L ndx 



A i-B 

n x D B d 3 Ldx + n x D%d 3 Ldx - 
J A 

B rb rb 

n AD x d±.Ldx + / n AD^d^Ldx — / n sD^d^Ldx 

1 JB JB 



n x D A d 3 Ldx 



r-A rB 

/ [ x D%d 3 L - x D a A d 3 L] ndx + / [d 2 L + x D a B d 3 L + A D^d 4 L] nd. 

Ja , J A 



[ A D p x d A L - B D^diL] r/da 



Since n is an arbitrary function, we can assume that n{x) = for all x £ [A, b] 
and so by the fundamental lemma of the calculus of variations, 



x D%d 3 L - x D a A d 3 L = 0, for all x £ [a, A}. 



Similarly, one proves the other two conditions: 

d 2 L + x D%d 3 L + A D%d 4 L = 0, for all x G [A, B], 

and 

A Dl d A L - B D^diL = 0, for all x G [B, b}. 

D 

3.2 The fractional isoperimetric problem 

The isoperimetric problem is one of the most ancient problems of the calculus 
of variations. For example, given a positive real number I, what is the shape 
of the closed curve C of length I which defines the maximal area? The most 
essential contribution towards its rigorous proof was given in 1841 and is due 
to Jacob Steiner (1796-1863). We state the fractional isoperimetric problem as 
follows. 

Given a functional J as in (O, which functions y minimize (or maximize) J, 
when subject to given boundary conditions 

y(a) = Va, y(b) = Vb, (10) 

and an integral constraint 

Hv) = f 9[y]dx - 1. (n) 

Here, similarly as before, we consider a function g of class C 1 , such that d^g has 
continuous right Riemann-Liouville fractional derivative of order a and d^g has 
continuous left Riemann-Liouville fractional derivative of order /3, and functions 
y G a-^b ■ A- function y G "E? that satisfies (jTUJ) and (fTTj) is called admissible. 

Definition 3.2. An admissible function y is an extremal for / in (jTTJ) if it 

satisfies the equation 

d 2 g[y](x) + x DZd 3 g[y](x) + a D^d 4 g[y](x) = 0, for all x G [a, b}. 



Observe that Definition 13.21 makes sense by Theorem 12.41 To solve the 
isoperimetric problem, the idea is to consider a new extended function. The 
exact formula for such extended function depends on y being or not an extremal 
for the integral functional I(y) (cf. Theorems 13.31 and 13. 4|) . 



Theorem 3.3. Let y be a local minimum for J given by {5[), subject to the 
conditions U0\) and ill\) . If y is not an extremal for the functional I , then 
there exists a constant A such that y satisfies 

d 2 F + x DZd 3 F + a D?d 4 F = (12) 

for all x G [a, b], where F = L + Xg. 



Proof. Let r\\,r\ 2 € Var(a, b) be two functions, t\ and e 2 two reals, and consider 
the new function of two parameters 



y = y + eirji + e 2 rj 2 . 



(13) 



The reason why we consider two parameters is because we can choose one of 
them as a function of the other in order to y satisfy the integral constraint. Let 



/(ei,e 2 ) 



g[y + eiTji + e 2 r] 2 ](x)dx -I. 



It follows by integration by parts that 

cb 



8I_ 

de 2 



(0,0) 



I (d 2 g m + d 3 g c a D% m + d i9 c x D?r, 2 ) dx 



(d 2 g + xDbd^g + aD%d 4 g) f] 2 dx. 



We have assumed that y is not an extremal for /, and therefore there exists a 
function r) 2 satisfying the condition 



dl_ 

de 2 



jLQ. 



(14) 



(0,0) 



Using (fi"4| and the fact that 1(0, 0) = 0, by the Implicit Function Theorem there 
exists a C 1 function e 2 (-), defined in some neighborhood of zero, such that 

J(ei,e a (ei))=0. 

Therefore, there exists a family of variations of type (|13[) which satisfy the 
integral constraint. 

We will now prove condition (|12[) . Similarly as before, we define a new 
function of two variables J(ti,e 2 ) = J(y). Since (0,0) is a local minimum of 
J, subject to the constraint 1(0, 0) = 0, and V/(0, 0) ^ 0, by the Lagrange 
Multiplier Rule (see, e.g., [32j p. 77]), there exists a constant A for which the 
following holds: 

V(J(0,0) + A/(0,0)) = 0. 

Simple calculations show that 



and 



dJ_ 
dl_ 



(0,0) 



(0,0) 



r b 

/ [d 2 L + x Dfd 3 L + a D%L] m dx 

J a 



[d 2 g + xD£d 3 g + a D^.d 4 g] rjidx. 



10 



In conclusion, it follows that 

cb 



f [d 2 L + x D£d 3 L + a D%d 4 L + X (d 2 g + x D?d 3 g + a D%d ig )] m dx = 0. 

J a 

By the arbitrariness of r\\ and the fundamental lemma of calculus of variations, 
one must have 

d 2 L + x D£d 3 L + a D^diL + X (d 2 g + x D£d 3 g + a D^g) = 0. 

This is equivalent to 

d 2 F + x D?d 3 F + JDgdtF = 0. 

□ 

Example 1. Let y{x) = E a (x a ), x £ [0,1], where E a is the Mittag-Lcmcr 
function: 

00 k 

E a (x) =^2 ^-, 7 , — , x£R,a>0. 



fe=0 



T(ak + 1) 



When a = 1, the Mittag-Leffler function is the exponencial function, Ei(x) 



The left Caputo fractional derivative of y is y (cf. [HI p. 98] ) , 

C a D a x y = y. 
Consider the following fractional variational problem: 

J(V)= f ^D a x yfdx -> extr, 

7i 1 

I{y) = j y^'D»ydx = J (y) 2 dx, 

k J/(0) = 1° and y(l) = E° a (l). 

The augmented function is 

F(x, y, ZD^y, c x D{y, A) = (% D«y) 2 + Xy^D a x y 

and the fractional Euler-Lagrange equation is 

d 2 F + x D«d 3 F + x D?d 4 F = 

i.e., 

x D«(2%D x 'y + Xy) = Q. 

A solution of this problem is A = — 2 and y = y. 

Observe that, as a — > 1, the variational problem (fT5|) becomes 

y' 2 dx — > extr, 



(15) 



11 



a = 0.1 




~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 

0.0 0.25 0.5 0.75 1.0 



a = 0.3 



a = 0.6 
a = 0.8 
a = 1 



Figure 1: Solutions of problem (|15 



yy'dx = -(e 2 - 1), 
/o z 

2/(0) = 1 and y(l) = e, 
and the Eulcr-Lagrange equation is 

d2 F _ ^Q 3 f = o & -±(2y' - 2y) = 0, 
ax ax 



(16) 



where F = y' 2 — 2yy'. Also, for a = 1, y(x) = e x , which is obviously a solution 
of the differential equation (fT6|) (cf. Figured]). 

We now study the case when y is an extremal of / (the so called abnormal 
case). 

Theorem 3.4. Let y be a local minimum of J ifS)J, subject to the conditions 
Iil0\) and ill]) . Then, there exist two constants Xq and X, with (Ao, A) ^ (0,0), 
such that 

d 2 K + x D%d 3 K + a D%K = 

where K = XqL + Xg. 

Proof. Following the proof of Theorem 13.31 (0,0) is an extremal of J subject 
to the constraint 1 = 0. Then, by the abnormal Lagrange multiplier rule (see, 
e.g., [351 p. 82]), there exist two reals Ao and A, not both zero, such that 

V(A o J(0,0) + A/(0,0)) = O. 



12 



Therefore, 



A 



dJ_ 



+ A 



(0,0) 






(0,0) 



The rest of proof is similar to the one of Theorem 



3.3 An extension 



□ 



We now present a solution for the isoperimetric problem for functionals of type 
|9]). Similarly, one has an integral constraint, but this time of the form 



i*(y) 



g{x,y{x), c a D«y{x),°Dly{x))dx = l. 



(17) 



Again, we need the concept of extremal for a functional of type (fT7|) . 

Definition 3.5. A function y is called extremal for I* given by (JT7J) if 

d 2 g[y](x) + x D B d 3 g[y}(x) + A D^d 4 g[y}(x) = 0, for all x G [A, B]. 

Theorem 3.6. If y is a local minimum of J* given by (0), when restricted to 
the conditions y(a) = y a , y{b) = yb and |_?7| j ; and if y is not an extremal for I* , 
then there exists a constant A such that 



d 2 F + x D B d 3 F + A D^d 4 F = for all x G [A, B] 
x D B d 3 F - x D a A d 3 F = for all x G [a, A] 

A D^d 4 F - B D l id A F = for all x G [B, b] 



(18) 



where F = L + Xg. 

Remark 3. In case [A, B] = [a, b], Theorem 13.61 is reduced to Theorem 13.31 

Proof. We consider a variation of form y — y + ei>7i + e 2 rj 2 , where 

T7i,??2 G c Var{a,b) and 771(A) = rji(B) = i l2 {A) = r) 2 (B) = 0. 
Define I* by the expression 

,-B 



/*(ei,e 2 ) 



g[y]dx - 1. 



Then, /*(0,0) =0 and 



dl* 
de 2 



(0,0) 



B r 



d 2 g m + d 3 g c a D a x n 2 + d i9 ^Dfo 



dx 



\ [ x D a B d 3 g - x D a A d 3 g] m dx + / [d 2 g + x D%d 3 g + A D^ x d 4 g] n 2 dx 

Ja , J A 



[ A D^d 4 g - B D^d 4 g] n 2 dx 



13 



(the last expression follows by integration by parts and some technical calcu- 
lations as presented in the proof of Theorem 13.11) . Let 772 be a function such 
that 

df * /o 



de 2 



(0,0) 



(its existence is guarantied since y is not an extremal for /*). Therefore, we can 
consider a subset of the family of functions {y + 61771 + £27/2 | (d, £2) & K 2 } that 
is admissible for the isoperimetric problem. Let J*(ei, £2) = J*(y)- Then, there 
exists a real A such that 



V(J*(0,0) + A/*(0,0)) = O. 



Similarly, one has 



dJ* 
~de\~ 



(19) 



(0,0) 



[ x D B d 3 L - x D a A d 3 L] Vl dx + / [d 2 L + x D a B dsL + A D%L] m dx 

J A 

b 

[ A D^d 4 L - B D^d 4 L] m dx, 



and 

dl* 

de 1 



(0,0) 



D B 8 3 g - x D a A 8 3 g) m dx + / [d 2 g + x D B d 3 g + A D^ x d 4 g\ mdx 

J A 

[ A D^d 4 g - B D^d 4 g] mdx. 



By©, 



dJ* 
~de\~ 



+ A 



(0,0) 



dl* 
~de\~ 



= 



(0,0) 



and ([T8l) follows from the arbitrariness of 771 . 



□ 



The following result generalizes Theorem 13.41 and is proved in a similar way. 

Theorem 3.7. If y is a local minimum of J* given by (0), subject to the bound- 
ary conditions \10\) and the integral constraint J_?7| ), then there exist two con- 
stants Ao and A, not both zero, such that 

d 2 K + x D B d 3 K + A D%d 4 K = for all x e [A, B] 
x D B d 3 K - x D a A d 3 K = for all x e [a, A] 

A D^d 4 K - B D^d 4 K = for all x E [B, b] 

where K = XqL + Xg. 
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3.4 Sufficient conditions of optimality 

We are now interested in finding sufficient conditions for J to attain local ex- 
tremes. Typically, some conditions of convexity over the Lagrangian are needed. 

Definition 3.8. We say that f(x,y,u,v) is convex in S C M 4 if d 2 f, d 3 f and 
dif exist and are continuous, and the condition 

f(x, V + Vi, u + Ui,v + Vi) - f(x, y, u, v) 

> d 2 f(x, y, u, v)y 1 + d 3 f(x, y, u, v)ui + d 4 f(x, y, u, v)v\ 

holds for every (x,y, u,v), (x,y + yi,u + u\,v + V\) 6 S. 

Theorem 3.9. Suppose that the function L(x,y,u,v) is convex in [a,b] x M. 3 . 
Then each solution yo of the fractional Euler-Lagrange equation $j) minimizes 

b 



J{y) = f L[y]dx 

J a 



when restricted to the boundary conditions y(a) — yo(a) and y(b) — yo(b). 

Proof. Let r\ £ "E^ be a function such that 77(a) = n(b) — 0. Then, using 
integration by parts, it follows that 

J(Vo + rj) - J(yo) = (L[y + rj\ - L[y ]) dx 

r b 

> 



I (d 2 L[y ] V + d 3 L[y ] c a D a x il + d i L[y^D^i) dx 

J a 

f [8 2 L + x D«d 3 L + a D^diL] [y ] ■ V dx = 

J a 



since yo is a solution of ((HJ. We just proved that J(y + rj) > J{yo)- D 

This procedure can be easily adapted for the isoperimetric problem. 

Theorem 3.10. Suppose that, for some constant X, the functions L(x,y,u,v) 
and Xg(x, y, u, v) are convex in [a, b] xl 3 . Let F = L + Xg . Then each solution 
yo of the fractional equation 

d 2 F + x D£d 3 F + a D%F = 

minimizes 

J{y) = I L[y]dx, 



under the constraints y{a) = yoifl) and y(b) = yo{b) and 

i{y) = / g\y]dx = i,leR. 

J a 
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Proof. Let us prove that yo minimizes 

F(y)= f (L[y] + Xg[y})dx. 

J a 

First, it is easy to prove that L(x, y, u, v) + Xg(x, y, u, v) is convex. Let rj £ "E^ 
be such that 77(a) = 77(6) = 0. Then, by Theorem [3^1 F(y + rf) > F(y ). In 
other words, if y g "^b is any function such that y(a) = yo(a) and y(b) = yo(b), 
then 

/>b rb />b />b 

/ L[y}dx+ \g[y}dx> L[y }dx+ Xg[y ]dx. 
If we restrict to the integral constraint, we obtain 



pb pb 

I L[y]dx + l> / L[y ]dx + l, 

J a J a 



and so 



no po 

I L[y]dx > / L[yo]dx, 

J a J a 

proving the desired result. □ 

Example 2. Recall Example [TJ Since L(x,y,u,v) = u 2 and \g(x,y,u,v) — 
—2yu are both convex, we conclude that y is actually a minimum for the frac- 
tional variational problem (1151) . 
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